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Abstract. We investigate the characteristics of FIR 
brightness fluctuations at 90 /im and 170 fim in the Lock- 
man Hole, which were surveyed with ISOPHOT aboard 
the Infrared Space Observatory (ISO). We first calculated 
the angular correlation function of each field and then 
its Fourier transform(the angular Power Spectral Density: 
PSD) over the spatial frequency range of / = 0.05 — 
1 arcmin -1 . The PSDs are found to be rather flat at 
low spatial frequencies(/ < 0.1 arcmin" 1 ), slowly decreas- 
ing toward higher frequencies. These spectra are unlike 
the power-law ones seen in the IR cirrus fluctuations, 
and are well explained by randomly distributed point 
sources. Furthermore, point-to-point comparison between 
90 yum and 170 /Ltm brightness shows a linear correlation 
between them, and the slope of the linear fit is much shal- 
lower than that expected from the IR cirrus color, and 
is consistent with the color of galaxies at low or moder- 
ate redshift(z < 1). We conclude that the brightness fluc- 
tuations in the Lockman Hole are not caused by the IR 
cirrus, but are most likely due to faint star-forming galax- 
ies. We also give the constraints on the galaxy number 
counts down to 35 mJy at 90 /im and 60 mJy at 170 /im, 
which indicate the existence of a strong evolution down 
to these fluxes in the counts. The galaxies responsible for 
the fluctuations also significantly contribute to the cosmic 
infrared background radiation. 
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1. Introduction 

In order to understand the history of galaxy formation and 
evolution, surveys at far-infrared(FIR) and submillimeter 
wavelengths are essentially important, since a large por- 
tion of star formation activity in the universe may be hid- 
den by dust, prohibiting optical and near-infrared studies 
due to the enormous extinction. The recent detection of 
the far-infrared Cosmic Infrared Background(CIB) radia- 
tion by COBE, which appears to have comparable bright- 
ness to the total intensity in dee p opt ical counts fro m the 
Hubble deep field (Fixen et al. tL998t Hauser et al. |1998; 



Lagache et al. 200C ), indicates that the infrared-bright 
galaxies are responsible for roughly half of the energy re- 
leased by nucleosynthesis. Various discrete source surveys 
are also being pursued from mid-infrared to submillimeter 
wavelengths. For example, source counts at 15 fim (Aus- 
sel et al. |1998|; Altieri et al. [199$ Elbaz et al. |1999|), at 



170 ^m (Kawara et al. |l998|(pape r I) ; Puget et al. |1999J), 
and at 850 /jm (Blain et al. |1999 ) have been reported. Us- 
ing both CIB and these new number count data, modeling 
of the cosmic star formation history has been attempted 



by many authors (Guiderdoni et al. 1998; Dwek et al. 1998 
Rowan-Robinson |1999[ Tan et al. |1999| ; Ishii et al. |1999[ ) 



All of them require a strong evolution in the star forma- 
tion rate as we look back to high rcdshift(more than 10 
times larger at z ~ 1). 

Although the FIR deep survey is now considered a 
key observing method for the exploration of the "opti- 
cally dark side" of th e sta r formation history of galax- 
ies (Guiderdoni et al. 1997 ), the detectivity of lm-class 



space FIR telescopes is likely to be limited by the noise 
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due to the fluctuation of the IR cirrus(Low et al. 1984), 
emission from the interstellar dust, even at high Galactic 



latitude(Hclou & Bcichman 1990). Hence a detailed study 
of the IR cirrus fluctuations is highly important, especially 
for the planning of the deep surveys intended with forth- 



and C_160(170/im). Each of the two fields is made up of 
4 sub-fields. 

Each sub-field map was produced from the edited raw 
data by the PHT Interactive Analysis(PIA; Gabriel et al. 



coming IR space telescopes such as IRIS (Murakami 1998) 
and SIRTF(Fanson et al. |1998| ). 

The spatial structure of the IR cirrus at 100 /im as 
measured by IRAS was extensively studied by Gautier 
et al. ( [1992D and Abergel et al.( |l996| ). From a Fourier 
analysis of the brightness distribution, Gautier et al. found 
that the Power Spectral Density(PSD) of the brightness 
fluctuation at 100 /im follows a power-law function of the 
spatial frequency with an index of about 3 below the spa- 
tial frequency corresponding to the IRAS beam size(about 
0.25 arcmin" 1 ). They also found that the PSD is propor- 
tional to Bq, where Bq is the mean b right ness of the IR 
cirrus. The ISOP HOT (Lemke et al. |l996| ) onboard the 
ISO(Kessler et al. 1996| ) is capable of observing the IR cir- 
rus at wavelengths longer than 100 /im, with a better spa- 
tial resolution than that of IRAS. Hcrbstmcicr et al. ( 1998 ) 
analyzed the spatial characteristics of four fields measured 
by ISOPHOT, and obtained similar power-law spectra for 
r elativ ely bright cirrus regions. Recently Lagache & Puget 
(2000) have made a power spectrum analysis of the 170 /im 
image of the Marano 1 field, and found a significant excess 
at / = 0.25 — 0.6arcmin _1 , which they attributed to the 
fluctuations due to unresolved extra-galactic sources. 

In this paper we investigate the characteristics of FIR 
brightness fluctuations in the Lockman Hole, a region with 
a uniquely low HI column density (Lockman et al. 1986| ), 
and thus minimal IR cirrus contribution to the fluctua- 
tions. Therefore the fluctuations in the Lockman Hole are 
likely to be dominated by faint, distant galaxies (Helou & 
Beichman 1990; Herbstmeier et al. 1998) and the fluctua- 
tion analysis of FIR images will provide unique informa- 
tion on the number counts of infrared galaxies even below 
the source-confusion limit, thus constraining the parame- 
ters characterizing the number count models. 

This paper is organized as follows: section 2 briefly 
describes the observations and data processing. Section 3 
explains the power spectrum analysis and examines the 
contribution by the IR cirrus fluctuations and also de- 
scribes a simulation of the images and the PSDs. Section 
4 describes the constraints on the galaxy number counts. 
The nature of the sources responsible for the fluctuations 
is discussed in section 5, and section 6 gives the conclu- 
sions. 



2. Observation and data processing 

The FIR survey of the Lockman Hole, which was executed 
as a part of Japan/UH cosmology program using the IS AS 
guaranteed time, is described in Taniguchi et al. (1994) 



1997) version 7.1 or 7.2, and is hereafter referred to as 



an AAP (Astronomical Analysis Processing of PIA) map. 
Each AAP map is either 58 x 58 pixels(23"/pixel) for a 
90 /im sub-field, or 31 x 31 pixels(46"/pixel) for a 170 /im 
sub-field. To correct the drift in the responsivity of the de- 
tectors, we applied the median filter smoothing (see Paper 
I for details) to the AAP data. 

Together with the final AAP maps we also produced 
the uncertainty maps, and found that the typical 3cr noise 
is as low as 0.012MJy/sr( 0.60mJy/pixel for a 46" x 46" 
pixel), indicating that the instrumental noise is negligible 
in the following results. 

As was done in Paper I, the observed fluxes as well as 
the brightness of the images are scaled based on the fluxes 
of the brightest source F10507+5723(UGC 06009) mea- 
sured with IRAS. We found that the flux calibration based 
on the FCS1 measurements underestimates the 90 /im flux 
of the IRAS source by a factor of 2.6, although it gives 
the mean brightness of the images consistent with the 
COBE/DIRBE brightness within 25 per cent. The 170 /im 
flux of the IRAS source is assumed to be 1133 mJy as de- 
scribed in Paper I. The flux calibration based on the FCS1 
measurements again underestimates the 170 /im flux by a 
factor of 1.5. For the 170 /im flux, we found that the dis- 
crepancy is mostly due to an underestimate of the effective 
solid angle of the ISOP HOT detector. Puget et al. fll999D 
and Lagache & Puget (2000) used an effective solid angle 
at 170 /im derived from Saturn footprint measurements, 
which is significantly larger than that used in the PIA. 
Here we assume theoretical PSFs of a telescope with a 60 
cm primary mirror and a 20 cm secondary mirror which 
agree well with those measured during the ISOPHOT cali- 
bration observations (Miiller 2000). With this assumption, 
the factors of 2.6(90 /im) and 1.5(170 /im) discrepancies re- 
duce to factors of 2.0(90 /im) and 1.2(170 /im). The origin 
of the residual discrepancies is still unknown. However, 
this is not problematic for the main results and conclu- 
sions of this paper as described in section 4 and 5 as long 
as the fluxes of F10507+5723 determined by the IRAS 
faint source survey are correct. 

Fig. [l] shows the images of LHEX and LHNW used 
for the fluctuation analysis, each of which is the largest 
square area extracted from the mosaiced map(see Fig. 
3 of Paper I) made up from four sub-field AAP maps. 
Each image is rebinned into 103 x 103 pixels 2 (C_90), 
53 x 53 pixels 2 (C.160 LHEX), or 51 x 51 pixels 2 (C.160 
LHNW). The plate scale is 23"/pixel for 90 /im maps, and 
46"/pixcl for 170 /im maps. 



and paper I. Two 44' x 44' fields named LHEX and LHNW 
were mapped with two filters: C_90(centered at 90 /im) 
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LHEX C_90 LHEX C_ 1 60 




LHNW C_90 LHNW C_160 



Fig. 1. The left column shows 90/zm im- 
ages of LHEX(top) and LHNW(bottom). 
Each image is 39'.5 x 39'.5 wide(23"/pixel). 
The right column shows 170/im 
images(46"/pixel). The LHEX image 
is 40'.6 x 40'.6 wide, while the LHNW one 
is 39'. 1 x 39'. 1 wide. The brightness shown 
in each image is offset from its median 
brightness. Each image is made of four 
22' x 22' sub-fields. Roughly the north is 
left and the west is top. 




3. The power spectra and the simulation 

For each image shown in Fig. [j] the 2-dimensional angu- 
lar correlation function C(x,y) was calculated from the 
brightness distribution B(xo,yo)' 

C(x, y) =< (B(x , y ) - B)(B(x +x,y + y)-B)> (1) 

where bracket < > represents the average over the whole 
area in each image shown in Fig. H, and B =< B >. 
Then, the coordinates were expressed in polar coordinates 
as C (r, 9) , and the Fourier transforms in various radial di- 
rections were calculated (i.e. 9 is fixed for each transform). 
One-dimensional power spectral density(PSD) P(f) is cal- 
culated from 

2 N ~ l 

P ^ = „ E C k co S (2nfr k ) , (2) 

fe=0 

where / is spatial frequency(l/(A^r p j x ) < / < l/(2r p i x )), 
r k = r p j X fc (k — 0, 1, • • ■ , N — 1), r p j X is the pixel size 
of the map, Nr p - lx is the largest angular size for which 
the angular correlation function is evaluated, and C k = 



C(r = rfe, 6) . Note that the units of the PSDs are the 
same as that of the angular correlation function: Jy 2 /sr 2 . 
Finally, the PSDs were averaged with respect to 9 and 
are shown in Fig. |2[ It is noteworthy that the fluctua- 
tions at high frequencies^ > 0.8 arcmin -1 for 90 /zm, and 
/ > 0.4 arcmin -1 for 170 /zm) are smoothed by the instru- 
mental beam and therefore the PSDs decrease appreciably. 
The small error bars of PSDs represent the standard de- 
viation among a set of the PSDs with different 9, showing 
that the PSDs are almost independent of 9. 

In order to check the contributions from bright sources, 
the PSDs are derived by masking circular regions with 
a 4xFWHM diameter around bright sources with fluxes 
above 5 max : S max = 250 mJy at 170 /zrr|] S max = 150 mJy 
at 90 /im. In the following the resultant PSDs are called 
"residual PSDs" , and are also shown in Fig. ^. Interest- 
ingly, the residual PSD for each image is more than half of 
the PSD of corresponding original image with almost the 
same spectral shape, indicating that there remains signifi- 

1 As for 170 fim images, smaller values of Smax significantly 
reduce the amount of residual images. 
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Fig. 2. Fluctuation power spectral densities(PSDs) of 
90/^m(top) and 170/im(bottom) images. Open circles rep- 
resent PSDs of LHEX and filled circles represent PSDs of 
LHNW. As well as the PSDs of the original images(Fig. |l]) 
shown by circles connected by thin-solid lines, the PSDs 
of the residual images ("residual PSDs"), where the pixels 
containing bright sources above 150 mJy(90/im), 250 mJy 
(170 jLim) are masked, are also shown by circles alone. The 
brightest source (IRAS F10597+5723) located in LHEX 
significantly contributes to the PSDs of the LHEX images. 
The lrj error bars, shown only for the residual PSDs, rep- 
resent the standard deviation in a set of the PSDs with 
different position angles in the sky. Thick line is an exam- 
ple of the IR cirrus PSD, which is an average spectrum of 
several IR cirrus in Ursa Major, with 100 /jm brightness 
of 2-3 MJy/sr. Dotted lines are the IR cirrus PSDs in 
the Lockman Hole, estimated by assuming that the cirrus 
PSDs are proportional to Bq, where B is mean brightness 
of the cirrus cloud. 



cant contribution from randomly distributed point sources 
with fluxes below S maDC . 

In Fig. H typical IR cirrus PSDs are also compared in 
order to check the contribution of the IR cirrus to the 



PSDs of the Lockman HoleQ We examined several IRAS 
100 (im maps of high-latitude clouds in Ursa Major(Z = 
145°, 6 = 40°), which are reproduced from the IRAS 
Sky Survey Atlas(ISSA) by reducing the brightness by 
a factor of 0.72, fo llowin g the COBE/DIRBE calibra- 
tion(Whcelock et al. 1994). The average brightness of the 
cirrus is 2-3 MJy/sr, and each map is 150' x 150' wide with 
l'.5/pixel. The derived PSDs show a power-law distribu- 
tion with an index of -1.5. Gautier et al.( |1992| ) noted that 
one-dimensional analysis of the IR cirrus yielded spectral 
indices near -2. These cirrus fluctuation spectra are much 
different from those obtained for the Lockman Hole images 
at 90 /im, which show rather flat spectra at lower spatial 
frequencies. The 170 /im spectra present a slope similar to 
the IR cirrus one, but this can be explained by the shape 
of the footprint power spectrum of ISOPHOT detectors. 
Moreover, the fluctuation powers are much larger than 
those estimated for the IR cirrus, which are also shown 
in Fig. H. We assume that the cirrus PSD is proportional 



to ^(Gautier et al. 1992), and taking the mean bright- 
ness Bo of the IR cirrus in the Lockman Hole as 0.33 
MJy/sr at 90 /im and 1.0 MJy/sr at 170 /im. These values 
are estimated from the atomic hydrogen column density of 
6 x 10 19 atoms citT 2 in LHEX(Jahoda et al. |199C| ) and the 
COBE/DIRBE data analysis in the Lockman Hole by La- 
gache et al. ( |1999| , in their Table 4). We can conclude that 
the IR cirrus contribution to the PSDs in the Lockman 
Hole is negligible over all spatial frequencies / > 0.05. 

In the following we interpret the residual PSDs in 
terms of unresolved point sources, probably galaxies. We 
neglect the spatial correlation between galaxies, thus as- 
suming that galaxies are randomly distributed in the im- 
ages. Then the residual PSD will be the product of the 
footprint power spectrum of the ISOPHOT detector and 
the fluctuation power due to the point sources, which is a 
white power spectrum given by: 



,dN 

Is' 



dS 



(3) 



where dN/dS is differential source counts. From the 
residual PSDs observed, we derive P S ourco = 13000 ± 
3000Jy 2 /sr at 90 /*m (S'max = 150 mJy) and P SOU rco = 
12000 ± 2000 Jy 2 /sr at 170 /im (S* max = 250 mJy), where 
the errors do not include systematic ones due to uncer- 
tainties in the flux calibration. Lagache & Puget ( 2000 ) 
reported the detection of P S ourcc = 7400 Jy 2 /sr at 170 fim 
for the Marano 1 field, after sources brighter than 100 
mJy are removed. Contribution from detected sources 
with fluxes between 100 mJy and 250 mJy is estimated to 
be about 7000 Jy 2 /sr. Thus the fluctuation power in the 



The formula for the IR cirrus PSD by Gautier et al. (1992) 



cannot be directly compared with the present work, because 
the power spectrum analysis presented here is one-dimensional 
one while the formula of Gautier et al. is based on the two- 
dimensional analysis. 
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Fig. 3. The residual PSDs of 90/im(top) 
170 /im(bottom) images(open circles and filled cir- 
cles, same as Fig. ^ are compared with the simulated 
PSDs based on various number counts models: dashed 
lines a re PS Ds of the simulated images by Guiderdoni 
et al. ( 1998 ) scenar io E, while the dotted lines are those 
by Takeuchi et al.( |l999| ) no-evolution. Thick gray lines 
are examples of the simulated images produced by simple 
double power-law number count models (see text for 
details). 



Lockman Hole due to the sources fainter than 100 mJy 
is ~ 5000 Jy 2 /sr, which is comparable to that observed in 
the Marano 1 field. 

A simulation was performed by making 90 /im and 
170 ^m images made up only by galaxies with fluxes be- 
tween S min and iSmax and calculating their PSDs. Here 
galaxies are treated as point sources with a PSF specific 
to the respective wavelength band of ISOPHOT. We used 
the image of the bright IRAS source (F10507+5723) seen 
in LHEX images as the PSF. The number of sources and 
their flux densities are controlled by the source counts. 
We examined the non-evolution count model by Takeuchi 
et al. (1999 ) and the scenario E by Guiderdoni et al. ( 1998] ) . 
We assume S m i n = lOmJy because the fluctuations due 
to galaxies fainter than 10 mJy are negligible in case of 
these models. The resulted PSDs are compared with the 



residual PSDs in Fig. |[ These simulated PSDs are not suf- 
ficient to explain the observed PSDs although the spectral 
shapes are quite similar. 

4. Constraints on the Number Counts 

In this section we will investigate the source counts below 
Smax, by fitting the simulated PSDs to the residual PSDs. 
In Paper I we presented the number density of sources 
brighter than 150 mJy at both filter bands. After Paper I 
was published, the data processing and the source extrac- 
tion technique have been improved, and now we obtained 
source number counts down to S = 70 mJy at 90 /im and 
S = 100 mJy at 170 /im, although the counts at such low 
fluxes may be overestimated due to the source confusion. 
We found that the cumulative counts above 150 mJy show 
quite a steep increase as the source flux decreases: (3 ~ 3 
where N(> S) oc . These results will be described 



in detail in Kawara et al. ( 2000 , paper III). If (3 remains 
equal to -3 down to S'min = 1 mJy, then the fluctuation 
power(P S ourcc) calculated by equation || exceeds the ob- 
served ones by an order of magnitude. Hence at a certain 
flux below S'max, we expect that the slope of the counts 
must flatten to {3 < 2 so that the predicted fluctuations do 
not exceed the observed ones. We thus considered simple 
double power-law count models and evaluated the simu- 
lated PSDs: in the flux Tcinge •S'max & ^crit 



N(> s) = N tt 



s 



and in flux range S C xjt ^ S ^ S n 



N(> S) = N crit 



Sent 



-Pi 



(4) 



(5) 



where iV crlt = iV max (S cri t /S'max) 00 , and /3 , ft are as- 
sumed to obey inequalities < ft < ft < 4. N max and 
its uncertainty are listed in Table |]. The uncertainty in 
JV max includes Poisson uncertainties based on the total 
number of sources with flux above S'max in all images, and 
systematic ones originating from the incompleteness due 
to the source confusion, which will be discussed in paper 
III. If Pi < 2 < /3q, the simulated PSD is not sensitive 
to Smin and is dominated by sources with fluxes around 
Scrit, as discussed in Lagache and Puget ( 2000| ). Various 
sets of parameters can be chosen so that the simulated 
PSDs fit to the residual PSDs. Examples of the simulated 
PSDs which fit well to the residual PSDs are also shown 
in Fig. |, for which /3 = 3.0, ft = 1.0 for both 90 ^m 
and 170 fim counts, and S^rit = 80 mJy for 90 ^m counts 
and Scrit = 150 mJy for 170 /im one. We found that the 
simulated PSDs can fit the residual PSDs over all spatial 
frequencies within ±22 per cent for the 90 /im PSDs and 
±17 per cent for the 170 /tm ones. These deviations of the 
residual PSDs from the simulated PSDs are also included 
as uncertainties in the residual PSDs. 
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Table 1. Number count parameters 



0.07 0.1 

S [Jy] 



Fig. 4. The shaded areas show the allowed regions for var- 
ious number count models, which are consistent with the 
fluctuation powers of the Lockman Hole images(Fig. ||). 
The bottom border is obtained from the simple num- 
ber count model with a single power-law index between 
Smin = 1 mJy and S max — 150 mJy(90/im) or 250 
mJy(170/xm). The upper border is determined from the 
allowed double power-law count models. See text for de- 
tailed information. The filled circles are observed num- 
ber counts at S max - The theoretical number count models 
by Guiderdoni et al.( 1998| ) (scenario E : solid lines) and 
the no evolution model by Takeuchi et al.( 1999| ) (dashed 
lines) are also shown. The open diamonds and triangles 
are observed and incompleteness-corrected source counts 
obtained by the FIRBACK Marano 1 survey(Puget et al. 
1999). The open circles with dash-dotted line show the 



IRAS 60 (im counts, in which the flux is scaled by the 
factor S v (90(j,m)/S u (60nm) = 2.1. 



We now derive the allowed regions in plots of the cu- 
mulative counts at fluxes below S max . The bottom border 
of the allowed regions is derived by examining the single 
power-law counts(S cr it = S m i n = 1 niJy)0 and we obtain 
Po=1.8 for 90 /jm counts and Pq=1.7 for 170 /zm ones. In 
case of double power-law counts we may choose a larger 
value of Po with a set of parameter pairs(S cr it, /3i), and 
among them S cr i t for /3j = is the largest. Thus we deter- 



band S'max 


max 


uncertainty 


[Jy] 


[sources/sr] 


in TVmax 


90 0.15 


1.4 x 10 5 


+ 16%, -47% 


170 0.25 


5.4 x 10 4 


+24%, -44% 



mine the upper borders of the allowed regions by connect- 
ing (S'crit, A'ci-it) points of numerous models with (3\ = 
and (3q = 2 — 4. The results are shown by the shaded area 
in Fig. |]. As for the 170 //m counts, the bottom border 
of the allowed region is close to the Scenario E model by 
Guiderdoni et al. 



(1998) 



and is consistent with the counts 
at 120 mJy and at 200 mJy ob tained from the FIRBACK 
Marano 1 survey (Puget et al. 1999). On the other hand, 



3 Lower value of S m 
suits: ,0o = 1.7 even if Smin 



does not appreciably change the re- 
0. 



the allowed region for 90 /jm counts are significantly above 
those of any currently existing models. 

The bottom border derived here cannot be applied 
below a certain flux Sl, which is the flux at the inter- 
section between an upper-border count model with the 
steepest slope (/3o = 4, /3i = 0), and the bottom bor- 
der: S'l = 30mJj/(90/xm) or 50 mJy(175^m). By defini- 
tion, the upper border model do not require sources below 
Scrit(> ^l)- O n t ne other hand, in case of the bottom bor- 
der model, contribution to the fluctuations from sources 
below S'l is still appreciable, about half of the total. Con- 
sidering the uncertainties in Af max and the residual PSDs 
as discussed above, we finally give Sl = 35 mJy for 90 /jm 
counts and Sl = 60 mJy for 170 /jm counts. 

5. Discussion 

We attribute the origin of the residual PSDs to the ran- 
dom distributions of faint sources which can no longer be 
identified as individual sources due to the source confu- 
sion. Another piece of evidence which supports this inter- 
pretation is the 170 /mi/90 /im brightness ratio of the FIR 
background emission. In Fig. |^ we show the point-to-point 
comparison between the 90 /im and the 170 /im sky bright- 
ness. Here we examine the residual images after masking 
pixels around the sources with S > S max . In order to 
match the spatial resolution between 90 /mi(76" FWHM) 
and 170/im(144" FWHM) images, one wavelength band 
image is convoluted with the other wavelength band beam 
profile. Although the scatter is large, the plot shows linear 
correlation with a slope of unity. This background color is 
quite different from that expected for the IR cirrus (about 
3.1), and can be interpreted as a typical FIR color of 
galaxies contributing to the fluctuation. The fluctuation 
color due to faint galaxies depends on their redshift, and 
their SEDs. The relations between the FIR flux ratio and 
the redshift for SEDs of the cirrus domi nated galaxy, the 
pure starburst ga laxy (E fstathiou et al. 200C ; Efstathiou 
& Siebenmorgen 2000 ), and the mixture of these two 
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Fig. 5. Point-to-point comparison between the 90 /im sky 
brightness and the 170 /im one in the Lockman Hole. 
The brightness is not absolute one, but is offset from 
the median brightness. In this plot the pixels around the 
bright sources^ > 150mJy for 90 /im, S > 250mJy for 
170 /im) are not shown. In order to avoid the error due 
to the difference of the beam size between the two bands, 
the 90 /im images are convoluted with the 170 /an beam 
profile(144"FWHM). The solid line showi ng the IR cirrus 
color vector is taken from Lagache et al. (199£), which is 
17. 5K gray-body with an emissivity proportional to A~ 2 . 



SEDs which represents a star-forming galaxy's SED like 
IRAS F10507+5723, are shown in Fig. 6. The flux ratio 
of 5^(170/Ltm)/S' 1/ (90/im) = 1 is attained at z ~ 0.7 for 
a pure starburst galaxy, while a cirrus dominated galaxy 
or a star-forming galaxy with a small contribution from 
starburst component must be local(z < 0.1). 

From Fig. 6 we can also justify the flux calibration 
based on the IRAS fluxes of F10507+5723. If we adopt 
the ISO 90 /im flux determined by PIA and the theoret- 
ical PSFs, then S' I ,(90/im)/5 I/ (60/im) = 1.0 which is the 
color of a pure-starburst galaxy, while the ISO flux ra- 
tio S' 1/ (170/im)/S , ^(90/im) = 1.6 is much larger than that 
expected for any nearby(z ~ 0) galaxies. 

The local cirrus dominated galaxies are not likely to 
be the dominant sources responsible for the fluctuations, 
since such galaxies are relatively bright at optical wave- 
length and hence the source number density in the Lock- 
man Hole at optical wavelength would exceed the ob- 
served one by more than an order of magnitude: as shown 
m Table @ we obtained source density of about 2 x 10 
sources/sr(60sources/deg 2 ) at S ~ 100 mJy. Using the 



(1997), the B magnitude of a 100 mJy source at 90 /im is 
estimated to be about 16 mag. Reported B-band source 
counts above B = 16 mag is only a few sources per square 



degrees (Kirschner et al. 1979; Ellis 1983). Hence the galax- 
ies responsible for the fluctuations must be heavily ex- 
tincted at optical wavelengths, which is a well-known fea- 
ture of starburst galaxies. In conclusion, the major source 
of the brightness fluctuations is most probably the star- 
forming galaxies located at z < 1 cxtinctcd at optical 
wavelengths. 

In the plot of the 90 /im counts in Fig. S, the 
IRAS counts are also plotted. These were originally 
the 60 /im counts, in which the flux is scaled by the 
factor S , ^(90/xm)/5'y(60/im)=2.1. This scaling factor is 
derived by the IRAS source counts at relativel y hig h 
fluxes (S„( 100 /im) > Uy) (Rowan-Robinson et al. 1986| ), 
and is consistent with the color of normal galaxies with a 
partial contribution from starburst. Thus, if such galaxies 
dominate the 60 /im counts down to 110 mJy, this scaling 
is valid. We used the complete differe ntial s ource counts 
down to H OmJy by Lonsdale et al. ( 1990| ) and Bertin 
et al.( 1997 ), and as shown in Fig. |J the scaled counts 
are in good agreement with the model counts given by 
Guiderdoni et al.(1998). Since the integrated counts down 
to 150 mJy is much larger than the model counts, the in- 
tegrated counts must show a steep rise between 150 mJy 
and 240 mJy, which, however, could not be observed in the 
IRAS 60 /im counts due to the source confusion (Hacking 
and Houck 1987). Further studies from space like IRIS and 
SIRTF will judge if this steep rise is real. 

Finally, we mention on the impact of the present 
work to the cosmic infrared background(CIB). Hauser 
et al.( 1998| ) reported an upper limit of 1.1 MJy/sr at 
100 /im and a positive detection of the CIB of 1.2 MJy/sr 
at 140 ^m and 1.1 MJy/sr at 240 (im. Lagache et al. fll999| ) 
reported lower values of the CIB: 0.72 ± 0.30 MJy/sr at 
140 /im and 0.91 ± 0.15 MJy/sr at 240 /im by further sub- 
tracting the dust emission associated with the diffuse ion- 



ized gas. Lagache et al.(2000) also reported the CIB at 



100 fim: 0.78 ± 0.21 MJy/sr. By summing the fluxes of 
all detected sources above 150 mJy, we obtained an in- 
tegrated brightness of the CIB of 0.031 MJy/sr at 90 /im 
and 0.050 MJy/sr at 170 /im. If we consider the constraints 
on the number counts for S > Sl shown in Fig. ^, then 
the integrated brightness would be 0.09 - 0.30 MJy/sr at 
90 /im for sources above 35 mJy and 0.053 - 0.15 MJy/sr 
at 170 /im for source above 60 mJy. Hence, 5-40 per cent 
of CIB can now be attributed to the integrated light of 
discrete sources above Si, which are responsible for the 
fluctuations. 

6. Conclusion 

We have studied the characteristics of FIR brightness fluc- 
tuations at 90 /im and 170 /im in the Lockman Hole, using 
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Fig. 6. FIR color vs redshift(z) relation for the cirrus dom- 
inated galaxy and the starburst galaxy. (Top) the ratio of 
5^(170 /urn) to S v (90 nm), (bottom) the ratio of 5,, (90 /im) 
to Su(60fim). T he mo del FIR SEDs are taken from Ef- 
stathiou et al. (2000) and Efstathiou & Siebenmorgen 
( 2000| ). The dash-dotted lines represent the color-;? rela- 
tion of an example of the mixed SED, cirrus : starburst = 
0.75 : 0.25 at 100 fan. 



the power spectrum analysis over the spatial frequency 
range of / = 0.05 — 1 arcmin -1 . The spectra of PSDs are 
found to be rather flat at low frequencies, which differs 
from that expected for IR cirrus fluctuations. We inter- 
pret the spectral characteristics of PSDs as those of ran- 
domly distributed point sources. The fluctuations in the 
Lockman Hole are not dominated by the IR cirrus, and 
are instead most likely due to star-forming galaxies. We 
next showed the constraints on the galaxy number counts 
down to 35 mJy at 90 /im and 60 mJy at 170 /im, assuming 
that only galaxies contribute to the observed fluctuation 
powers. This analysis indicates the existence of steep rise 
in the integrated counts down to these fluxes and espe- 
cially at 90 /im, the source density is much larger than 



that expected from the currently available number count 
models. 

We also found a linear correlation between sky bright- 
nesses at 90 /xm and 170 /im, even if pixels around the 
bright sources are masked from the images. The slope of 
the linear fit is much lower than the color of the IR cirrus, 
indicating that the point sources are star-forming galaxies 
at low or moderate redshift(;z < 1). These galaxies with 
fluxes above 35 mJy or 60 mJy also significantly contribute 



to the CIB recently reported by Hauser et al.(1998) and 
Lagache et al.( |2000[ ). 
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